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Abstract. 

We investigate the time evolution of a model system of interacting particles, moving in 
a (i-dimensional torus. The microscopic dynamics are first order in time with velocities set 
equal to the negative gradient of a potential energy term \& plus independent Brownian 
motions: \& is the sum of pair potentials, V(r) +j d J('yr), the second term has the form of 
a Kac potential with inverse range 7. Using diffusive hydro dynamical scaling (spatial scale 
7 _1 , temporal scale 7 -2 ) we obtain, in the limit 7 | 0, a diffusive type integro-differential 
equation describing the time evolution of the macroscopic density profile. 
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1. Introduction. 



The transition from the microscopic dynamics of interacting particles to hydrodynam- 
ical type equations describing the coarse grained evolution of macroscopic variables, such 
as the diffusion equation for the density, is a basic problem of non-equilibrium statistical 
mechanics. While far from resolved for systems with realistic interactions there has been 
much progress recently on this problem for model systems. Like in real systems, the tran- 
sition from microscopic to macroscopic evolutions in these models is based on a separation 
between microscopic and macroscopic scales. Setting e equal to the ratio of microscopic 
to macroscopic spatial scale and then looking at macroscopic times which are of order e~ a 
microscopic time units, a = 1 for Euler (non-dissipative) and a = 2 for diffusive evolutions, 
we expect to obtain the macroscopic equations in the hydro dynamical scaling limit (HSL) 
e I 0. We refer to the books of De Masi and Presutti, [3], and Spohn, [24], for a general 
background on this subject (see also the review article by Lebowitz, Presutti, and Spohn, 
[16]). 

To actually prove this HSL, one needs to show that during macroscopic evolutions the 
microscopic particle system can be well described, on the microscopic scale, by a local ver- 
sion of the equilibrium measure which is stationary under the dynamics. These measures 
depend on quantities conserved by the microscopic dynamics, such as the particle density, 
which then evolve on the slower hydrodynamic time scale according to the hydrodynamic 
equations. This requires good mixing or chaotic properties of the dynamics (as well as of 
all the relevant equilibrium states). This is particularly so for the case of diffusive scaling 
where longer times are involved. It is for this reason that the only model systems of inter- 
acting particles for which the HSL has been established in the diffusive limit are systems 
with stochastic dynamics. Thus the HSL for Ginzburg-Landau models was established first 
by Guo, Papanicolaou and Varadhan, [12], by applying entropy techniques. These tech- 
niques were further developed by Rezakhanlou, [21], to cover the case when the invariant 
measure is not a product measure and phase transitions may occur. These methods can be 
applied also to lattice gas models that satisfy the so called "gradient condition", [24]. For 
lattice gases this condition is however not natural and the only known examples are when 
the invariant measure is a product measure or the spatial dimension is one, [13], [24]. Very 
recently the diffusive HSL for non gradient lattice gases has been proved by Varadhan and 
Yau, [26]. For systems of particles in the continuum the gradient condition is more natu- 
ral, while a common technical problem in these models is the control of the local number 
of particles: the conservation law cannot prevent locally very high densities. The only 
continuum models treated with the entropy techniques quoted above are one dimensional 
systems of Brownian particles interacting via positive superstable short range potentials 
considered by Varadhan, [25], and Ornstein-Uhlenbeck interacting processes studied by 
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011a and Varadhan, [19]. We should also mention here that the diffusive limit can be 
proven for a Hamiltonian system of non-interacting particles moving among a fixed array 
of convex hard scatterers: the Sinai billiard system with finite horizon in d = 2, [1], [2], 
[17], [18]. 

In 1991 Yau, [27], proposed a new method for proving the HSL of interacting particle 
systems of gradient type, looking at the relative entropy and its rate of change w.r.t. local 
Gibbs states. This method can be applied also to continuum systems in higher dimension, 
e.g. in the derivation of the Euler equations from a Hamiltonian system with weak noise 
considered by 011a, Varadhan, and Yau, [20]. 

In the present paper we extend the work of Varadhan to Brownian particles with positive 
superstable short range potentials in all dimensions. In addition we also permit long range 
pair interactions of the Kac type in which the range parameter 7 -1 goes to infinity as the 
macro to micro spatial scale e~ l . This extends previous work for such systems on a lattice, 
[9]. 

To be more precise, we consider a system of N particles which evolve in time according 
to the non-inertial Brownian dynamics 

dr- f)^> 

_i = _|_ (ri ,..., rK ) +W , W p.!) 

where Wi(r) is a stochastic Langevin force with Gaussian statistics having covariance 
(P/2)8ijS(T — r')l, 1 the unit <i-dimensional tensor. The parameter f3 is the inverse tem- 
perature of the canonical ensemble, \i ~ exp[— (3^], which is the stationary measure for 
the evolution. The potential energy \1/ is a sum of pair potentials, 

• • • , r N ) = [ V ( r v) + -r dj (Vij)] (1-2) 

where ty, = 7\ — rj and the r\, i = 1, . . . , N, are confined to a <i-dimensional torus of 
length L. We take V(r) to have a finite range R with R < 7L: 7 -1 is the range of the Kac 
potential which will be taken to be large compared to the inter-particle spacing L/N 1 ^. 
Systems with interaction of form (1.2), with V(r) = 0, J(r) > 0, and different types of 
dynamics, have been investigated numerically and analytically by Klein and coworkers as 
model of glassy dynamics, [11], [14], [15]. 

We observe that due of to the prefactor r ) d+1 appearing in the force term due to the 
Kac potential, the dynamics defined by (1.1) is a weak perturbation of the one defined for 
J = 0, i.e. without long range interactions. Thus we may expect that for small 7's the 
system reaches local equilibrium w.r.t. the short range potential on spatial scales smaller 
than 7 _1 at times of order 7 -2 . The effect of the long range interaction on such states 
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will then appear only in determining the macroscopic equation for the relevant parameters 
describing the local equilibrium. 

In fact we shall take as our initial distribution something close to the local equilibrium 
distribution relative to the short range potential V with a density which varies on the scale 
of L ~ N 1 ^ ~ 7 _1 and consider macroscopic times of order r/7 -2 . The HSL will then 
correspond to letting 7 J, 0. We will prove that in that limit the density profile on the 
macroscopic scales x and t will satisfy the following non-local integro-differential equation 
of the diffusive type: 

^(t, x) = V • x))Vp(t, x) + a(p(t, x)) j^dy VJ(x - y)p(t, y)| (1.3) 

where the integral is over the <i-dimensional unit torus T d and a(p) = f3p is the mobility 
of a system of interacting Brownian particles, which, due to the fact that the system is 
gradient, does not depend on the interactions, [24]. The diffusion coefficient D(p) is given 
explicitly in terms of the Helmholtz free energy density a(f3, p) associated to the "reference 
system" interacting only with the short range potential V, in such a way that the following 
"Einstein Relation" holds ([24]): 

D{p) = a{p) Tp =a{p)^ (1.4) 

where A is the chemical potential of the reference system at density p. As in the lattice 
case, [9], eq. (1.3) can be rewritten in terms of the gradient flux associated to the classical 
local mean field free energy functional and the density dependent mobility <r(p): 

^(t, x ) = V-L(p)V 5 -f)(t,x) (1.5) 



dt v ' ' \ y " bp 
where 

F{p)= / dxa(p,p(x))+ - dx dy J (x - y) p(x) p(y) (1.6) 
J T d Z J T d J T d 

Our proof is based on Yau's method quoted above. The main restriction of this method 
is that the derivation of the HSL is valid only as long as the macroscopic equation has a 
smooth classical solution. Consequently, unlike the lattice case, we can no longer guarantee 
existence of global solutions. In fact, even if the initial datum is smooth and lies in the 
one phase region (for the reference system), we cannot guarantee that the time evolution 
will not develop singularities or create regions of high density where the reference system 
undergoes a phase transition and the diffusion coefficient D(p) vanishes. 
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The outline of the rest of the paper is as follows. In Section 2 we give a precise description 
of our system and present the results. In Section 3 we prove the HSL by computing the 
relative entropy and its rate of change w.r.t. the local equilibrium states of the reference 
system. To do this we need a local ergodic theorem whose proof is sketched in Section 4 
and large deviation estimates for the local Gibbs states which are the content of Section 
5. A local existence theorem of classical solutions for the macroscopic equation is quite 
standard, a sketch of the proof is given at the end of Section 3. 



2. Notation and results. 

In this section we state our problem in a precise mathematical form using from the 
beginning the rescaled space and time variables, xi = 7^, and t = 7 2 r. We also ab- 
sorb (3/ 2 into the Brownian motion term which remains invariant under this rescaling of 
space and time. In these units we consider a system of N interacting Brownian motions 
x(t) = {x\{t), . . . , £jv(£)} with state space T d , the <i-dimensional unit torus, satisfying the 
following equations (i = 1, . . . , N): 



dxi = —(3 



dt + V2dwi (2.1) 



where {wi, . . . ,wn} are independent Brownian motions on T d , the parameter f3 > is the 
inverse temperature, J G C 2 (T d ) and Vir) e C 1 (IR d ) is a positive function of |r|, with 
compact support and such that V(0) > 0. The latter implies that V is superstable. In eq. 
(2.1), V7(7 _1 (ii — xj)) and VJ(xi — Xj) are the gradients of the functions V(-) and J(-) 
w.r.t. their arguments, evaluated at the points r y~ 1 (xi — Xj) and (x{ — Xj) respectively. 
We shall further assume that the number of particles iV depends on the scaling parameter 
7 G (0, 1] in such a way that N^ d / 1 as 7 [ (typically N = [y~ d ]). 
The process t — > x(t) is a diffusion on T dN with generator 

L 7 = 4°) + U 7 (2.2) 

where 

4 0) = E A * " T-^W^- 1 ^ - x 3 )) ■ V, (2.3) 

i i=lj 

and 
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In (2.3) and (2.4) A; (Vj) denotes the Laplacian (gradient) w.r.t. the i-th particle com- 
ponent of x G T dN . Note that the diffusion L 7 is reversible w.r.t. 



l^-y(dx) = — exp 



cte (2.5) 



where Z 7 is the normalization factor making p^ a probability measure on T dN . 

If the initial distribution of the diffusion has a density w.r.t. p,^ then the density 
at any later time, f-y(t,x), satisfies the forward Fokker-Planck equation 

where L* is the adjoint of L 1 w.r.t. // 7 . 

To state our result we need to introduce some thermodynamic quantities relative to 
the reference system, i.e. the system of particles interacting only via the short range 
(superstable) potential V. For any regular domain A of IR d we define the grand canonical 
partition function 



Z A (P, A) = e-' A l -Jrr / dri ' ' • drN exp " ? £ V ^ - r ^ 



^V=0 AN 



(2.7) 



where A G iR is the chemical potential. The pressure is defined by the limit 

p()M)= lim rf ^TA7 lo g Z A(/3,A) (2.8) 

which exists and defines a convex and continuous function of /5 and A, see [22] and [23]. 

Setting the inverse temperature equal to some fixed value j3 > (which we will some- 
times omit) there exists, for the reference system, a non empty open set U C iR such that 
for any X E U there is a unique (infinite volume) Gibbs state. This is a point process on 
iR d , invariant and ergodic w.r.t. space translations, satisfying the DLR equations relative 
to the potential V, see e.g. [6], [22]. The pressure is a smooth function of A G U and 
the average density of particles p, as a function of the chemical potential, is given by the 
smooth 1-1 map A i— > p(X) = d\p((3, A) of U onto W = d\p(/3,U) C iR + . To make more 
symmetric the correspondence between the parameters A and p we introduce the Helmholtz 
free energy a{f3, p) as the Legendre transform of the pressure: 

a(J3,p)= sup{Ap-p(/3,A)} (2.9) 
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and we recover the chemical potential as a function of the density by the smooth 1-1 map 
p i— > X(p) = d p a(/3, p) of W onto U. 

We consider the nonlinear non-local integro-differential equation (1.3) that we rewrite 
below in a more concise form: 



where '*' denotes convolution on T d and recall that cr(p) = (3p. In the sequel we will use 
the capital letter P to denote the pressure as a function of the density. Then P'(p) = pA'(p) 
so that the diffusion coefficient in (2.10) is D(p) = (3P'(p) (see (1.4)). 

In the one phase region the pressure is a smooth, strictly increasing function of the 
density, so that D(p) is smooth and strictly positive for any p G W. Then the following 
theorem holds, whose proof is sketched at the end of the next section. 

Theorem 2.1 There exist locally classical solutions of (2.10) that lie inside the one phase 
region W. 

We fix such a solution p(£, x) , < t < T (T > 0) . We may assume that there is a compact 
set K w C W such that p(t, x) G K w for any (£, x) G [0, T] x T d and diat(K w , 2R+\W) > 2<Ji 
for some <5i > 0. Clearly X(t, x) = d p a((3, p(t, x)) lies in the compact set K u = d p a((3, K w ) 
and dist(K u , M\U) > 25 2 for some 5 2 > 0. 

We introduce the local Gibbs state associated to the above macroscopic evolution p(t, x) 
as the probability measure on T dN which is absolutely continuous w.r.t. p 1 (dx) with 
density 



where C 7 (t) is the normalization constant making / 7 a probability density. 
Our main result is 

Theorem 2.2 Let / 7 be the solution of the Fokker-Plank equation (2.6) with an initial 
distribution such that 



dp 
di 



(t, x) = V • {D{p)Vp + (T(p)VJ * p}(t, x) 



(2.10) 




(2.11) 




(2.12) 



Then, for any (p G C°°(T d ), any 5 > 0, and any t G [0,T], 




(2.13) 
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where 



N i y2 { P( x i) ~ / dx(p(x)p(t,x) 



> 5 



Notation: From now on we will write /(£, •) = /(£) for functions on [0, T] x T d or [0, T] x 

j-dN 

We will prove Theorem 2.2 by using the relative entropy method introduced by Yau, 
[27]. We recall the basic entropy estimate: if \i,v are two probability measures on the 
same measurable space, then for any F € L x (dv), 



J dpF < H(p\v) + log J dv exp[F] 



where H(p\v) is the relative entropy of \x w.r.t. v and, if fx <^ u, 

dp 



H(fi\u) = J dp log 
For any t G [0, T] define the functional 



dv 



H 1 {t)= 1 d fp 7 (dx)f^x)\ogfy^- 



(2.14) 



Note that H 1 {t) is 7 d times the relative entropy of / 7 (£)<i/i 7 w.r.t. / 7 (t)(i / u 7 and that the 
argument of the limit in the l.h.s. of (2.12) is exactly if 7 (0). In the next section we will 
prove: 

Theorem 2.3 Under the same hypothesis of Theorem 2.2, for any t e [0, T], 

hm# 7 (t) = (2.15) 



The hydrodynamic limit (2.13) follows as a corollary of Theorem 2.3. To see this we 
note that it follows from the large deviation principle (LDP) for the local Gibbs states 
(2.11), see Section 5, that there is a c(8, (p) > such that 

iE^ (t) [l A t J < exp [ - c(5, (p)N] 

where IE^[-] denotes the expectation w.r.t. the measure jdii 1 and lr is the characteristic 
function of the set V. On the other hand, from the basic entropy estimate the following 
inequality holds (see e.g. [27]): 

mMt) log2 +7 -^ W 

'- J " log (l + ^ <'»[!,, J" 1 ) 

so that, for some C > 0, iE^ (t) [l A t ] < C(N~ 1 + (A^ 7 d )" 1 if 7 (t)) -> as 7 | 0. 
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3. Proof of Theorems 2.3 and 2.1. 



Because of the hypothesis (2.12) on the initial distribution, we only need a good estimate 
on the time derivative of H^(t). By Lemma 3.1 of [20], the following bound holds: 

Recalling (2.2) and that L 7 is reversible w.r.t. fi^(dx), (3.1) gives 
dH. 



dt 



(3.2) 



17 < l d f ^{d^f^^f^x)- 1 ^ - §^Ut,x) 

+ 7 d J n 7 (dx) / 7 (t , x) / 7 (t , x) ~ 1 17* / 7 (t, x) 
where U 1 is defined in (2.4). By an explicit computation 

f 7 (t, x)- 1 (l^ - ; 7 (*, x) = p {^1 VA i 2 ^ ^) + AA ^ ^) - A ^ x <) 

-7" 1 5^ /3Vy( 7 - 1 (x l -x J ))-VA(t,^)|+iE^W ^/3A(t,^) 

(3.3) 

where A denotes the time derivative of A. Integrating by parts one computes the action of 
the adjoint operator U* and gets 



-7" 1 /3V^(7 _1 (^ - x k )) -VJ( Xl - Xj )\ 



(3.4) 



In both (3.3) and (3.4) there is a term of the following type: 

K v (x) = ^ 7 " 1 VF( 7 - 1 (^ - x k )) ■ Vif^xi) 

for some smooth function (p on T d (actually <p = /3 2 A(t, •) and <p = /3 2 J(- - Xj) in (3.3) 
and (3.4) respectively). Since is an odd function, 

K v {x) = \Y, 7" 1 VF( 7 - 1 (^ - x k )) ■ (V^(xi) - V^(x fc )) 

(3.5) 



where is the £-th component ofi 6 T rf and D^ v = d 2 / (dx^dx v ). Since W has compact 
support, we can estimate the reminder using Taylor expansion: 

\R v {x)\ < r 7 (p) |Vy|(7 _1 (^ - x k )) 

with r-y(<p) — > as 7 | 0. 

Inserting (3.3) and (3.4) into (3.2) and using (3.5) we get 



dH. 



dt 



1 <lE f ^ t) [/3^ 1 (x,X(t))] +iE^ (t) ^20\(t,Xi) 



+ £ 7 (£) 



(3.6) 



where 



^(x, \(t)) = 7 d E [® VA i 2 ^ ^) + AA ^' ^) " ^) } 







+ ^ E E ^ (V^^T" 1 ^ - *;))0 C „A(f , xOCt" 1 ^ - xj))" 



+ l 2d E {PV\(t, Xi) ■ VJ{ Xi - X j) + AJ( Xl - Xj)} 



+^ 2d E E E^ w ) ? ^" 1 ^-^))^^-^)(7- 1 (^-^))" 



(3.7) 



i^j k:k^i £,77 

while £ 7 (£) satisfies the bound 

|e 7 (f)| < 7 d K(/? 2 A(t)) + A^r 7 (/3 2 J))iE^ | W^V* - x 3 )) 

To obtain the behavior of e 7 (t) when 7 | we use the following lemma which is proved in 
Section 4: 

Lemma 3.1 If W is a continuous function on M d with compact support, there is Cw > 
such that, for any 7 e (0, 1] and any t G (0, T], 



7 d E^ - v*-^)) 



By applying the lemma with W = \ W\ we conclude that 

lim sup |e 7 (t)| = 
Ti°te[o,T] 



(3.8) 



10 



Moreover, by the LDP for the local Gibbs state, see Section 5, for any t E [0, T], 



lim IE^ {t) 

7.10 



y X(t,Xi) = / dx A(£, x)p(t, x) 
a J Jr d 



(3.9) 



Now we want to write $ 7 Qe, X(t)) in terms of local empirical quantities. Let O be the 
space of particle configurations on IR d , i.e. w G O is a subset of IR d which is locally 
finite (see Section 5 for more details). Given x E T d , for any x E T dN we construct a 
configuration uj Xj1 E O by setting 



to. 



:,7 = Ui = 7 - x) : |Xi -x\< 1/4} 



(since there is no risk of confusion, to simplify notation we omit the explicit dependence 
on x of oo xa ). Clearly u Xyl is well defined in every compact set inside the cube of IR d of 
side 1/(27) an d centered in the origin. So, if F(u) is a local function on Q, F(u xa ) is well 
defined for any 7 small enough. 

Let us introduce the cubes D n = {q E M d : \q^\ < n, £ = 1, . . . , d}, n E IN. For any 
local function F we denote by F n its spatial average over the cube D n , i.e. 



F n (oj) = 



D 



dr F(r r u>) 



(3.10) 



n\ J D 



where r r is the space translation by r (r r r' = r + r'). 

Let x be any non negative function on IR d with compact support and total integral 1. 
We define the following local functions on O: 



(3.11) 



G*»=f E x(q)(VV)t(q- q ')(q- q r, 77 = 1 d 



(3.12) 



q,q Euj 



and let R n (uj), G^{uj) be their averages over D n . Observe that R(u>) is a natural version 
of local density for the configuration u, while G^ v (uj) is the local quantity appearing in the 
virial theorem (see below, before (3.19)). 

Lemma 3.2 Let cp, ip be smooth functions on T d . Then 



lim sup lim sup sup IE^^ 

n^oo 7 J.O te[0,T] 



7 d V"^(xi) -/ dxip(x)R n (oj xn ) 
Jr d 



(3.13) 
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lim sup lim sup sup 

n^oo 7 J.O te[0,T] 



i d E § (v^Ht- 1 ^ - x i ))( 7 - 1 (x 4 - Xj ))M 



.1, ; 



lim sup lim sup sup 

n^oo 7J.O te[0,T] 



- / dx<p(x)Gto(u; x „) =0 



- dx dyip(x)R n (uj xn )ip(x -y)R n (u yn ) 
J T d J T d 



lim sup lim sup sup IE^^' 

n^oo 7 J.O te[0,T] 



(3.14) 



(3.15) 



7 2d E E f(v^(7- 1 (^-^))(7- 1 (^-^)r 



X 



ifj(xi-Xj)- dx dyGi v (u Xyl )t(j{x -y)R n (u ya ) 
J T d J T d 



= 

(3.16) 



Proof. We have 



where, for z <E T d , 



/ <ix^(a:) J R n (w X)7 ) = 7 d E^, 7 (^i) 



V?n,7<»=7 d / dx<p(x)—— I dqxil 1 (z-x) + q) 

JT d \ u n\ J D n 

Then the expectation in the l.h.s. of (3.13) can be bounded by \\ip — </? n , 7 ||oo that vanishes 
as 7 I for any n G IN because of the smoothness assumptions on tp. In an analogous way 
one can estimate the expectation in the l.h.s. of (3.14) by 



\<P - ¥V 7 ||oo7 



d Jf]fi (t) 



Eliwi(7- 1 (^-^))|7- 1 (^-^)| 



and (3.14) follows from Lemma 3.1. 
Let us consider now (3.15). We have 

dx dyLp(x)R n (u Xjl )ij(x-y)R n (u y ^)=7 2d y^Lp n}1 (x i )'ifj n)1 (x i -x j )+Err 1 
J T d J T d *rr. 

where -i/> n , 7 is defined as ip nn and 

\Err t \ < ||v?n, 7 ||oo max{|^ nj7 (x - z) - ip n „(y - z)\ : x, y, z E T d , \x-y\< >y(r + n)} 
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with r such that the support of \ is contained in the closed ball of radius r. Since 

|V>n, 7 (z - z) - ipn^iy - z)\ < ||V^|| 00 7(r + n) + 2||V> - Vy 7 lloo, #m vanishes as 7 | 0. 
Since \ip nn (xi)ip nn (xi - xj) - ip(xi)ip(xi - xj) \ < \\<p\\oo\\tp ~ VV7IU + IMIoo||<P - ¥V 7 lloo, 
(3.15) follows. In the same manner we compute 



dx dyGi v (uj X;1 )ij(x-y)R n (ujy tl ) 
j T d j T d 



with 

\Err 2 \ < max{\ip n ^(x - z) - ip n ^{y - z)\ : x,y,z E T d , \x - y\ < *y(r + n)} 
x 7 d ]T IVFKt" 1 ^ -*,-)) |7 _1 (*i - Xj )\ 

Then (3.16) follows from Lemma 3.1. □ 

Collecting together (3.6), (3.8), (3.9) and applying Lemma 3.2 to IE^ [$ 7 (x, \{t))] , 
we obtain 



lim sup lim sup sup 

n^oo 7J.0 te[0,T] 



dH, 



dt 



1 _|_ jEf-zit) 



P I dx (^ x (t, x) - A(£, x)p(£, x)) 



< (3.17) 



with 



A(t, x) - /3| VA| 2 (t, x) - AA(£, x) R n (cv x , 7 ) - ^ D iv X(t, x)G^{u xn ) 
/3VA(t, x) • (V J * # n (u;., 7 )) (x).R n (u; X)7 ) - (A J * R n {u. n )) (x)R u (uj Xj1 ) 



where, for any <p G C(T d ), 



(3.18) 



(<^*i? n (o;. ;7 ))(x) = / dx(p(x-y)Rn(u y „) 

JT d 



Now we want to substitute the spatial average G^{u xn ) with a function of the empirical 
density R n {ixs xn ). More precisely we would like to replace it by the average of G^ v w.r.t. 
the Gibbs state with density equal to R n {uj xn ). 
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To do this we need to introduce some cutoffs. Let K be a compact set such that 



K w cKcW, dist(K,m+\W)>8 1 , dist(K w ,M+\K)>8 1 

(recaff that K w is the compact set inside the one phase region W where the solution p(t, x) 
lies and that dist(K w , IR + \ W) > 25i) and define the local function u n {uj) = l.K(Rn(u)). 
We denote also by (frk the cutoff at the level k G IR+, i.e. (j>k{s) = s if \s\ < k, 4>k{s) = 
sign (s) k otherwise. Finally let G^ v (p), p G W, be the average of G^ v (u) w.r.t. the unique 
Gibbs measure with density p. By the virial theorem, see e.g. [25] , 

G^{p) = ((3P{p)-p)kv 

where P(p) is the pressure as a function of the density p introduced just after (2.10). 
For any measurable function m :T d — > IR + we define the functional 



(3.19) 



fi(t, x, m) = X(t,x) -/3|VA| 2 (t, x) Jm(x) — /3AA(t, x)P(m(x)) 

- /?VA(£, x) • (VJ * m)(x)m(x) - /3(AJ * m)(x)P(m(x)) 
Observe now that, since P'(p) = pA'(p), by integration by parts, 

/ dxP(p(t,x))AX(t,x) = - dxp(t,x)\'(p(t,x))Vp(t,x)-V\(t,x) 
Jr d Jr d 

= - dx p{t,x)\V\{t,x)\ 2 

and, analogously, 

/ dxP(p(t,x))(AJ* p(t))(x) = - dxp(t,x)V\(t,x)-(VJ*p(t))(x) 
Jr d Jr d 

so that, for any t G [0, T], 

/ dx A(t, x)p{t, x) = I dx 0(t, x, p(t)) 
JT d Jr d 

Then we can replace A(t, x)p(t, x) by 0(t, x, p(t)) in (3.17). 
We decompose now, for any k > 0, 

4 

x) - 0(t, x, = °A x ) 
p=i 
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with 

CLi(t, x) = [Cl(t, x, R n (u. n )) - tt(t, x, p(t))]u n (u xn ) 

Q 2 (t, x) = [*W(t, x) - n(t, x, Rn(w.„))] 

n 3 (t, x) = [*W(t, x) - n(t, x, p(t))} (1 - u n (u; x „)) 

Cl 4 (t,x)= [^Jt,x)-^(t,x)] 

where ^ ] (t,x) is defined as Vx(t t x) in (3.18) with Gfp replaced by (fa o G^), 
In Section 5 we will prove that there is 5q > such that 



(3.20) 



lim sup lim sup lim sup sup \ JE^ ^ 

k^oo n^oo 7 J,0 te[0,T] 



dx (3Q p (t, x) 



(3.21) 



On the other hand, the local ergodic theorem, see Section 4, implies that 

rT 

7J.0 Jo UT d 



lim sup lim sup lim sup / dsIE^^ 

k^oo n^oo j[0 JO 



dx /9|0 2 (s, x) 



= 



(3.22) 



From (3.17), (3.21) and (3.22) we get, for any t G [0,T], 



H y (t) + rfs|iE^ (s) J dxf3VL 1 {s,x) - 25 ~ 1 H 1 (s)^ < 0(71,7) (3.23) 



with 



lim sup lim sup o(n, 7) = 

n^oo 7J.O 



Now, from the basic entropy estimate, for any S > and any s G [0, T] , 



dx (3Qi(s, x) 



5^- d [ dxpSlifax) 



so that, from (3.23), for any t G [0, T], 



H 7 (t) - <T V / ds log exp 



Ms) 



5-f~ d / dx(Xl!(s,x) 

Iq-d 



-(2d Q - 1 + d~ 1 ) [ dsH y (s) <o(n, 7 ) 



By applying the Gronwall Lemma to the last inequality we get 



H 7 (t) < e (2S o 1+s W 0(71,7) +<TV / ^ logiE^ (s) exp 



5<y d I dx fifths, x) 

Iq-d 



(3.24) 
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In Section 5 we will prove that 



lim sup lim sup 5 log iE^ 7 ^ exp 

n^oo 7J.O 



-5-f~ d [ dxf3VL 1 {s,x) 



<8- L e 5 (s,\) (3.25) 



where 



©^(s, A) = sup { I dx 



6(3 [0(s, x, m) — Q(s,x, p(s))] l^(m(a;)) — Ip(\(s, x),m(x)) 



m :T d — > iR + integrable 



and (recall definitions (2.8) and (2.9)) 

i/?(A, m) = A) + a(/3, m) - Am) 

From (3.24) and (3.25), for any t G [0,T], 

-t 



limsup# 7 (t) < e (2<5 < 



o 1 ^- 1 ^ [ dsS-'Qsis.X) 
Jo 



(3.26) 
(3.27) 

(3.28) 



We conclude the proof of Theorem 2.3 from (3.28) by showing that, for 5 small enough, 
@d(s. A) = for any s G [0, T]. We first note that, for A G U, m 1— > Ip(X,m) is strictly 
convex on if, non negative, and equal to iff m = d\p((3, A), so that 

— / dec Ip(\(s, x), m{x)) < 0, = iff m(er) = p(s,x) 
J 7 > 

On the other hand the functional 

Q s {m) = / dx (3 [0(s, cc, m) — 0(s, cc, p(s))] l^(m(cc)) 

is bounded on the class of functions considered in (3.26) and equal to for m = p(s). 
Then, for 5 small enough, 65(5, A) = provided that 



SQs_ 
5m 



(p(s)) = 



(observe that p(s) is away from IR + \ K because dist(K w , IR + \ K) > 8\). By an explicit 
computation, 

/T 1 H (p(s)) = X(s) - p(\ VX\ 2 (s) + AX(s)P'(p(s)) - VA(s) • V J * p(s) 



VJ * (p(s)VA(s)) - P'(p(s))AJ * p(s) - AJ * P(p(s))) 
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But, recalling that P'{p) = pX'(p) = (5- 1 D(p), 

| VA| 2 + AXP'(p) = X'(p) (Vp • VA + pAX) = /?- 1 A'(p)V • (D(p)Vp) 



and 



VA ■ VJ * p + VJ * (pVA) + (AJ * p)P'(p) + (AJ * P(p)) 

= A'(p) (Vp • V J * p + pAJ *p) = A'(p)V • (pVJ * p) 

So that, for any (s, x) E [0, T] x T d , 



^ (p(s)) iA'(p(.s..r)) 



^-V-{/J(p)Vp + /3pVJ*p} 



(s,x) = 



since p(s,x) satisfies (2.10) and a(p) = (3p. □ 

We conclude the section with the proof of Theorem 2.1. 

Proof of Theorem 2.1 (sketch). Let p E C 2 (T d ) be such that p (T d ) C W. By continuity 
and compactness we can find three compact intervals Ii = [aj,6j], i = 1,2,3, such that 
Po(1~ d ) Q h C h C I3 C W, a 3 < a2 < ai < b\ < 62 < We construct two functions 
J},aG C 1 (iR) with the following properties: D(u) = D(u) and a(u) = a(u) for -u E J 2 , 
c _1 < L>(-u) < c for some c > 1 and for any u E IR, supp(a) C I 3 . Then we consider the 
Cauchy problem 

^(t,x) = V- {D{p)Vp + a(p)V J * p}(t, x) (3.29) 

with initial datum po. Arguing exactly as in Theorem 4.1 and Remark 4.1 of [10], we 
know that there exists a (unique) classical solution p(t, x) of the Cauchy problem above 
(moreover it lies in the region J 3 at any time). Clearly we can find T > such that 
p(t, x) E h for any t E [0, T] and x E T d . From the choice of D and a it follows that 
{p(t, x); (t, x) E [0,T] x T d } is also a (local) classical solution of the original equation 
(2.10). □ 



4. Local ergodicity and entropy bounds. 



We start this section by proving (3.22). Since D^ V X and D^ V J * R n (u. n ) are bounded 
functions on T d , it is sufficient to prove that 



lim sup lim sup lim sup / dtlE^^ / dx | (4>k G^ v ) n (uj Xil ) 

k^oo n^oo 7 J,0 JO \_JT d 



G^{R n {u x ^))\u n {u x ^) 



(4.1) 



= 



The main step in proving (4.1) is a local ergodic theorem for the measure / 7 (t)<ip~ 
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Theorem 4.1 For any local, bounded and continuous function F, 

lim sup lim sup / dtlE^^ I dx \F n (u; X:1 ) - F(R n (uj X:1 ))\u n (uj X:1 ) 

n^oo 7J.0 JO \_JT d 



(4.2) 



where F{p), p G W, is the average of F w.r.t. the (unique) Gibbs measure with density 
P- 

Proof. We introduce the translation invariant density 



fj(x) = ^ ds dxfj(s,T x x) 
so that (4.2) is equivalent to proving that 



(4.3) 



lim sup lim sup JE^< 



= 



(4.4) 



Let Hp G -CV be such that supp(F) C -D nF (-D n as in definition (3.10)) and let n = n + np- 
Through the mapping x 1— > cuo,7, for any 7 small enough it is well defined the projection 



U fl :T dN ^Q\ D 



= ^0,7 In 

1 ±y r 



We have to characterize the family Tn of the (weak) limit points of \y 1 = H n (f^dpj); 7 G 
(0, 1]}. First of all we observe that, by translation invariance, for any finite region A of 

IE"-* [N A ] = iV 7 d |A| < |A| (4.5) 

This proves that the family \y 1 = n^(/ 7 <i/i 7 ); 7 G (0, 1]} is tight. 

Let jtt^n jf j be the canonical Gibbs measure on the cube D n with boundary conditions 
Co G O and number of particles m. We prove that any v G T n can be written as 



u(du) = J v(duj, dm) ^ n {du) 



(4.6) 



where v(duj,dm) is a measure supported on {m < \D n \}. Since any limit point satisfies 
an inequality like (4.5), we have only to prove (4.6) for some v. By a straightforward 
extension to higher dimensions of the analogous argument in Varadhan, [25, Lemma 7.5], 
we can reduce the problem to the estimate of a certain Dirichlet form. More precisely (4.6) 
follows if 



f lV-f I 2 

sup 7 d / df^ V 1 V < +00 

76(0,1] J • J 7 

The bound (4.7) is consequence of the following lemma: 
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(4.7) 



Lemma 4.2 Let 



f^(t,x)=- [ ds [ dxfj(t, t x x) 

t JO JT d 



and define, for any density f, 



a 7 (/) = 7 d J d^ 



/ 



Then there is C > such that, for any 7 G (0, 1] and any t G (0, T] 

<ry(f>r(t)) < t 



(4.8) 



(4.9) 



Proof. With an abuse of notation, denote by H(fj(t)\l) the relative entropy of / 7 (£)<i/i 7 
w.r.t. dfj,^. Observing that L 7 1 = we get, after some standard computations, 



^H^m) = j t Jd lh f 7 (t)]ogf 7 (t) = |^ 7 L;/ 7 (t)log/ 7 (t) 



(4.10) 



Since for any x G T d r x f 7 (t) = / 7 (£, r x -) solves the same Fokker-Planck equation, recalling 
definition (4.3), from (4.10) we get 

H(TMt)\l)-H(T x h(0)\l) 



dx 



- \ I ds j dx-f d a 1 (r x f 1 (s)) 
t Jo Jr d 

- J d// 7 7 d l3VJ{ Xl - X j) ■ Vi/ 7 (t) 



(4.11) 



Now, since cr 7 (-) is a convex functional, [4], 

<r 7 (/ 7 (t)) < - / cis / dx <r^(T x fi(s)) 
1 Jo Jr d 

On the other hand, by Cauchy-Schwartz inequality, 
- J d^ 7 d ^/?VJ(x, - xj) • Vi/ 7 (t) 



(4.12) 



< 



2 



(4.13) 
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with C\ = /3||VJ||oo- Also, since dfi y is t x -invariant, 

H(r x f 7 (0)\l) = ff(/ 7 (0)|l) = 7" d ^ 7 (0) - logC 7 (0) 

+ / ^ 7 / 7 (0,x)^/3A(0,x,) < ~f- d H 7 (0) + 2/3||A(0, -)||ooiV < C 27 ~ d 

J i 

(4.14) 

for some C 2 > (in the first bound we used Jensen's inequality, in the second one the 
assumption (2.12) on the initial distribution). Collecting together (4.11), (4.12), (4.13) 
and (4.14), recalling also that the relative entropy is a positive function, we get 



— f T d < -I'^iUt)) + Ci7" V ff 7(/r(*)) 
so that, for any 7 G (0, 1], 



a 7 (/ 7 (t))<C lV / a 7 (/ 7 (t)) + ^ (4.15) 

But (4.15) implies that <t 7 (/ 7 (£)) < C/t with C the positive solution of x = CiVTx + C<z- 
The lemma is proven. □ 

Now we conclude the proof of Theorem 4.1. Using (4.6), the l.h.s. of (4.4) can be 
bounded by 



limsup sup IE^ 



\F n {u)-F{R n {uj))\l K {R n {uj)) 



where Gi is the class of Gibbs states with density p < 1. The characteristic function 1 K 
reduces the problem to computing the above limit in the one phase region. The limit is 
then zero by the law of large numbers for the unique Gibbs state of given density p G K. 

Finally, the limit (4.1) follows easily from Theorem 4.1. In fact (4.2) with F = oG^ 
implies that the l.h.s. of (4.1) can be bounded by 

lim T sup {U k Td^(p) - G^(p)\} (4.16) 

fc^oo peK 

But, for any p G K, 

<f>k^G^(p) ~ G^{p) = JE^ [(cf> k o G^)(fj) - G^(u)] 

where p p is the (unique) Gibbs state with chemical potential A = d p a(/3,p). We observe 
now that (ftkoG^ 71 — > G^ v pointwise as k — > 00. Moreover \(j)k°G^ v \ < \G^ V \ < cNg for some 
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c > and some finite subset B of IR d . Recalling that, by superstability, JE^ P [Ng(u)] < oo, 
the limit (4.16) is zero by the Dominated Convergence Theorem. 

We conclude the section by proving Lemma 3.1: 

Proof of Lemma 3.1. Since V is positive and superstable, there is a constant Cw such 
that, for any x G T dN , 

Wfr-^xt - x 3 )) <C W J2 V^ix, - x 3 )) 

The previous inequality is a straightforward extension to higher dimensions of the anal- 
ogous one derived in the proof of Lemma 4.2 of [25], so we omit the details. Then it 
is enough to prove the lemma for W = V. From the basic entropy inequality, recalling 
definition (2.5) of the reference measure /U 7 , 



2 7 d f 
< — — log J dx exp 







^5>(7-Vi-^)) 



and, by Jensen inequality, 

2 7 d 



P 



log J dx exp -^Y^V^-^Xi-Xj)) <7 d jdxY v (r<~\ x i- x j)) < W V 



Then we are left with an estimate of the relative entropy H(f y (t)\l). Recalling (4.10) and 
that <r 7 (-) is a positive functional, we can bound 

tf(/ 7 (t)|l)<tf(/ 7 (0)|l)- 7 d fds [dfj, y Y^pVJixi-x^-Vif^s) (4.17) 

Jo J ijtj 

Since ^(dx) is translation invariant, recalling (4.8) and using (4.13), we have 

-l d J ds j d^ ^/3VJ(x, - Xj ) • Vi/ 7 (s) < C 1 t 7 - d v /a 7 (/ 7 (t)) (4.18) 

The r.h.s. of (4.18) can be bounded using (4.9). Then, recalling (4.14), from (4.17) and 
(4.18) we finally get, for some C > 0, 



2 7 c 



ff(/ 7 (f)|l)<C 



(4.19) 



The lemma is proved. □ 
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5. Large deviation estimates and removal of the cutoffs. 

Part of the large deviation estimates of this section are contained in the theory developed 
in [7], [8] and [20]. We will sometimes refer to these papers for proofs and details. 

Let us start with some elementary facts in the theory of point processes. A configuration 
of particles in IR d can be represented by a locally finite subset u of IR d . Sometimes it can 
be useful to look at w as a Radon point measure on IR d via the map oj \— > ^2 qeto S(- — q). 
We denote by Q the set of all such configurations. O can be made into a Polish space under 
the vague topology tq, defined as the smallest topology making continuous the mappings 
lo i— > joo{dq)g(q) = J2 q eu>9(l) f° r an y 9 '■ ~^ ^ which is continuous and compactly 
supported. The natural cr-algebra T on O is the one generated by the counting variables 
Nb ■ oo — > card(u; fl B) for B any Borel subset of IR d . It can be proven that T is the Borel 
cr-algebra relative to the topology tq . 

A point process on M d is a probability measure Q on (fi,jF). We denote by M. the set 
of all point processes Q with finite expected number of particles IE^[Nb] in any bounded 
Borel set B C IR d . M. can be equipped with the topology t w of weak convergence based 
on the topology tq. However it is useful to introduce a finer topology on A4, called the 
topology T£ of local convergence. Let £ be the class of measurable functions F on O 
that are local and tame, i.e. for any such F there are a bounded set B and a constant 
c> such that F(u) = F{oo n B) and |F(u;)| < c(l + N B (u)). Then r £ is defined as the 
weak* topology on M. relative to £, i.e. the smallest one making continuous the mappings 
Q l— *■ Q(F) = IE [F] for any F G C Note that tc is strictly finer than t w , as follows 
observing that the mappings Q h- > Q(Nb) are T£-continuous for any bounded Borel set B. 

Let M T be the set of all the stationary point processes Q in M, i.e. those Q such 
that Q = Q o t~ 1 for any x G M d . M T is r^-closed and is assumed equipped with the 
induced topology. For any Q G M. there is a Radon measure vq(cIx) on iR d such that 
Q{N B ) = vq{B) for any Borel set B C iR d . If Q G M T then u Q (dx) = p(Q)dx for some 
positive number called the intensity of Q. 

Let P be the Poisson point process on M d with intensity p{P) = 1 (i.e. P is such that for 
any collection of disjoint bounded subsets Bi, . . . , B n , the counting variables , . . . , Ns n 
are independent and Poisson distributed with parameters \B\\, . . . , |S n |). We introduce 
the entropy density h(Q\P) of Q G M. T w.r.t. P as follows. For any B C iR d we denote by 
the projection on S, and let D n , n G iAT, be as in definition (3.10). Denote by Ti n (Q\P) 
the relative entropy of Q o w.r.t. P o . It is easy to see that n i— > 7Y n (Q|P) is a 
super-additive functional, so we can define 

h(Q\P) = lim \D n \- l H n {Q\P) = sup |D n |- 1 W n (Q|P) (5.1) 

n^oo n 

Let now 1/ be a positive and superstable finite range potential like the one introduced 
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in Section 2. The associated Hamiltonian in D n with free boundary conditions is 

Hn{") = \ E ^(9"?'), ^^0 (5.2) 

For each n e IN and Q G -M r let 

be the expected energy per volume in D n . By our assumptions on the potential, $ n is well 
defined and positive. Moreover, see [8, Thm. 1], the limit &(Q) = lim^oo $> n (Q) exists 
and satisfies 

m) = lQ(U) XQeM T , 2 (53) 

^ +oo otherwise 

where 

M r ,2 = {Q e Mr : Q{N%) < +oo for any bounded B C M d } 

and 

UW = l E x(q)V(q-q') (5.4) 

In (5.4) x is an y non negative function on IR d with compact support and total integral 1 
(it is easy to verify that Q(U) does not depend on x if Q is a stationary point process). 
Finally define, for any A G IR and j3 > 0, 

Kp, x (Q) = mQ) + KQ\p) - PMQ) (5-5) 

The following proposition establishes the basic properties of the functionals introduced 
above (see [7], [8] and [20]). 

Proposition 5.1 The functionals Q,Kp\ : A4 T — > [0, +00] are lower semicontinuous 
relative to T£ (and then also relative to any coarser topology on M. T , e.g. t w ). Moreover 
Kp t \ has Tc-compact level sets. 

We can formulate now the LDP for Gibbsian point processes. The form we need here 
is strictly contained in the results of the papers quoted above: 
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Theorem 5.2 Let F G C and F n be as in definition (3.10). Then, for any A G IR and 
(3 > 0, there exists the limit 

lim lig- 1 logiE p exp [(3XN Dn - (3H n + \D n \F n ] = sup {Q(F) - Kp \(Q)} (5.6) 
and the r.h.s. of (5.6) is finite. 

Remark. In [8] the principle is formulated in terms of functionals of the stationary empirical 
field R n ,uj obtained by replacing the configuration u with the periodic continuation u/ n ) of 
its restriction to L> n , but the same result holds also for the spatial average F n by standard 
arguments that we do not describe here, see e.g. [7]. 

An immediate consequence of the above principle is the existence of the pressure p(f3, A) 
and a variational formula for it. In fact, taking F = in (5.6) and recalling (2.7) and 
(2.8), we get 

p(P,X) = - min f5- 1 K^ x {Q) (5.7) 
QeM T 

and, comparing with (2.9), we recover the Helmholtz free energy as 

a(P,p)= min {$(Q) + P^HQlP)} (5.8) 

QEM T :p(Q)=p 

where we used the fact that the r.h.s. of (5.8) is a convex function of p (which is true since 
$ and h(-\P) are affine functionals on M T ). 

Now we discuss the LDP for local Gibbs states of the type defined in (2.11). We forget 
the dependence on t G [0, T] and we deal with a generic smooth map x i— > X(x) of T d into 
the one phase region U such that the corresponding density x i— > p{x) G W satisfies 



/ dxp(x) = l (5.9) 



We denote by f-y(x) the density w.r.t. p^dx) of the associated local Gibbs state (observe 
that the functions x \— > A(t, x) introduced in Section 2 satisfy the conditions above for any 
t G [0, T]). The following theorem is contained in Section 5 of [20]. 

Theorem 5.3 Let A G £ and (p G C(T d ). Then 



lim 7 d log /dx exp V /3A(x,) - J V V(7 1 (x i -x J ))+7 d / rfxv?(x)A(w X) 
= SUP { f rd dX [v( X )Qx( A ) ~ K (3,\(x)(Qx)]; {Qx} c M T ■. j^dxp{Q x ) = 



(5.10) 



24 



Taking A = in (5.10) and recalling (5.7), we get 
lim7 d log Jdxex.p ^ f3\(xi) - - ^ Vfr -1 ^ - x 3 )) 



(5.11) 



= snp {J Td dx [~ K 0Mx)(Q*)]'> {^} : J^ d dx P(Qx) = l| = J^dxpp(p,X(x)) 

where in the last equality we used the fact that the supremum of the integrand is reached 
at p(Q x ) = p(^) satisfying (5.9). 

Theorem 5.3 is not sufficient for our purposes since to prove (3.25) we need also an 
upper bound for the Laplace asymptotics of non local functionals. The following theorem 
gives the required estimate. 

Theorem 5.4 Let A,F,Ge C with F bounded and let <p E C(T d ) and E C(T d x T d ). 
Then: 

lim sup 7 d log IE f ~< exp7~ d / dx<p(x)A(u; x „) + dx dyip(x,y)F((ju x ^)G(u y ^) 

7 J.O \_JT d JT d JT d 

< sup I J^dx J^dy [(p(x)Q x (A) + ^(x, y)Q x (F)Q y (G) 
-Kp Mx) (Q x )-pp(p,\(x))];{Q x }dM T : [ dx p(Q x ) = 1 

JT d 



(5.12) 



Proof. For any configuration of particles 00 on T d define 

T{u,) = Y,PK*)-^ E V( 1 -\z-z'))+ 1 - d j dx<p{x)A{u; x „) 



z,z 6w 
z^z' 



(5.13) 



+ 7 d I dx I dyip(x,y)F(u Xi y)G(uy,y) 
J T d J T d 



Setting 



Z C N = Jdx exp[T(x)] 
and making use of (5.11) it is enough to prove that 

limsup7 d logZ^ < supi I dx I dy [cp(x)Q x (A) + i/;(x,y)Q x {F)Q y (G) 

7 10 I JT d JT d 

-Kp Mx) (Q x )] - {Q x } c M T : J^dxp(Q x ) = lj 



(5.14) 
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First of all we observe that for any partition of T d into cubes of side 2e there are 
uniform approximations of A, cp, and that are constant on this partition. Let T e (-), Z^ e 
be defined as T(-), with A, and ^ replaced by their approximations. Let -D^ be a 
cube (as in definition (3.10)) such that the functions A, F and G depend only on u H Di 
and, for some c > 0, A, G < c(l + Nd 1o )- Then one easily bounds 

|T(x)-T e (x)| <0!(e) 7 - d / dx(l + N D (u x „)) < 1 (e)\D £o \(l + N) (5.15) 

with Oi(e) — > as e J, (uniformly in 7), so that 

lim sup 7 d log Z C N < lim sup 7 d log Z C N e + O t (e)\D to \ 

7 10 7-L0 

On the other hand, also the error made in replacing the r.h.s. of (5.14) with its e- 
approximation (i.e. the one defined with the piecewise constant approximations of A, <p, 
and ip) is easily bounded by 2 (e)|-D^ H | with 2 (e) — > as e j 0. Thus, with no loss of 
generality we can prove (5.14) assuming \,<p, and ip constant on some cubic partition of 
T d 

Let P 1 be the Poisson point process on T d with intensity r )~ d and define, for any \x G IR, 



d ^ ^-dn ef in 



Z„( 7 ) = ^ exp[T(u;) + = £ ^— — ^ 

Since 

1 , e-^iV^ 
lim — log — — = 

then (recall that iV 7 d / 1 as 7 j 0) 

lim sup 7 d log < inf limsup { 7 d log Z M ( 7 ) — A*} (5.16) 
740 ' M^/zf 7 ^o 

From (5.16) we get (5.14) if 

lim sup 7 d log ^(7) < sup dx dy {(p(x)Q x (A) + tp(x, y)Q x (F)Q y (G) 

740 {Q*}cM T JTd J T d (5.17) 

- K(3,\(x)(Qx) +MQx)} 

By redefining X(x) as X(x) + fx, it is enough to prove (5.17) for [i = 0. Divide T d into 
disjoint boxes B a of center a and side (2 7 £) and let _D^ be the cube as in (5.15). Set 

B aM = {x E B a : dist(x, T d \ B a ) > 2 7 £ } (5.18) 
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and define 



T aa ,{oo) = {2 1 ty 



z,z €.u>nB a 
z^z' 



+ 7 d ip((T,a') / dx dy F(u x „)G(u}y„) 



'tr,Z - cr' ,e 



(5.19) 

Since the potential is positive, we obtain an upper bound on T(u) by neglecting the 
interaction between different boxes. Then 



T{UJ) < YjTaa'i^ + Kiuj) 



(5.20) 



where 1Z takes into account the errors due to the integration on the smaller cubes B a ^ in 
(5.19). Because of the choice of £q we easily bound, for some C > 0, 

^<C 7 " d J dx{l + N Dlo {u xn )) 

U(B CT \B CT ,, ) 

We decompose \J(B a \ B a ^ ) into a union of disjoint cubes of side (27^0) • Calling x a 
the centers of these cubes, we have 



n „ 

n<Cj2 dr(l + N DlQ (r r+7 -i 4a o;o l7 )) 



with n ~ (7^) d (£/£o) d 1 so that, by Jensen inequality, 



lim sup lim sup 7 d log JE P ^ exp [7£] < lim sup lim sup 7 d n log iE p exp [(7(1 + N De )] = 

(5.21) 

where we used the independence and stationarity of P 1 and that IE P exp[aNs] < 00 for 
any a > and any bounded set B C iR d . From (5.20) and (5.21) we get 



lim sup 7 d log -Zq (7) < lim sup lim sup 7 d log IE P "< exp 

7J.O £->oo 7 ].0 



(5.22) 



Now define the following local function on fi x O: 

S aa , (wi, c^ 2 ) = pXi^RM - pufa) + ip^AM + ^(a, </)F(u;i)G(u;2) 
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where R and U are denned in (3.11) and (5.4) respectively. Expanding variables in the 
r.h.s. of (5.22) and neglecting errors that can be proved to be of order 0(1/ 1) uniformly 
in 7 G (0, 1] as done in (5.21), one easily obtains (introducing a new parameter e = 2£^) 



limsup7 logZ (7) < limsuplimsup X(i, e) 

7J.0 l^oo ej.0 



(5.23) 



with 



X(£,e) = e d \ ^| _1 log IE®° P ° exp 



e d ^2\D £ \(S <ja ') £ (oj a ,oj a ') 



0,0' 



(5.24) 



where ® a Pa is the product measure on x • • • x Q of e d independent Poisson 

processes P a on JR d and, analogously to (3.10), for any local function M on O x O, we 
defined 

Mg(uJi, uj 2 ) = 777-m dr dr' M(r r u;i, Tt-'U^) 

1-^7 I J D e J D e 

To analyze the limit in the r.h.s. of (5.23) we need first to introduce some cutoffs. Given 
k + ,k- G IN, let Xk + k_(£,e) be defined as X(£,e) in (5.24) with S aa > replaced by 



k + if S aa ' > k + 
S h „J - = { S aa , if -k_ < S aa , < k 
-k- if S arT i < —k- 



+ 



Then, by definition of relative entropy, 

X k+k _(£,e) < sup ilE® \e 2d J] (S^-)^,^)] - e d \D t \~ 1 H(Q\ ® a (P a o 7r^ + _)) j 

(5.25) 

where the supremum is taken over point processes Q on .Dg + ^ x • • • x D e+ j>. The parameter 
£ is chosen so large that S a(7 ' is D^-measurable. The variational formula for the relative 
entropy gives also the explicit form of the measure Q where the supremum in the r.h.s. of 
(5.25) is achieved, see e.g. [5, Prop. 1.4.2]: Q <C ® a (P a o 7r^ _) and 



dQ 



= M e exp 



(5.26) 



(H e the normalization constant). For any x G T d let a(x) be such that x G B a ^, and 
denote by Q a , Q a a> the a-th, {a, a'}-th marginals of Q. From (5.26) it is easy to prove 
that, for any x, y G T d , i^t/, 

kmP>v & v(^Qa(x)a(y)(dui,di02), J Q(dui)Q a{x ) (dui \^{a{x)}")Qa{y) (^2 \^>{a(y)}-)^ = 



28 



where D var (-, •) denotes the variation distance between measures. Moreover 
H(Q\ ® CT (-Per o n~ D ]J) > J2 H (Q°\ P ° ° *d)J 

a 

Since A, cp, ip are piecewise constant, as e is small enough, for any x,y E T d , the function 
S a (x)a(y) is independent on the partition {B a } and equal to 

S xy {u u u 2 ) = f3X{x)R{u 1 ) - /3U (wi) + + ^(x, y)F(a; 1 )G'(a; 2 ) 

Thus the limit as e I of the r.h.s. of (5.25) is easily bounded and we get 

limsup* fc+fc _ (*,e)< sup / dx f dy{Q' x ®Q' y {S k x f-)-\D,\- 1 H{Q' x \Po^-) )\ 

(5.27) 

where the supremum is taken over all the collections {Q' x ;x G T d } of point processes on 
D £+ i. Now the proof follows in a standard way, see e.g. [7] , [20] . We extend Q' x to a point 
process Q x on M d by taking independent copies on all disjoint cubes translated of D i+ g. 
Then we obtain from it a stationary process Q x ^ by setting 

= \D t ^T x I drr r Q" x 

J D e+e 

From convexity of the relative entropy and the independence properties of P one easily 
proves that h(Q x \P) < \D ^i^ 1 H (Q X \P o 7r^* _). Moreover, for any bounded local 

function MonOxfi, \Q x e) <g> Q y e) (M) -Q' x ® Q' y (M e )\ -> as £ -> oo. Then from (5.27) 
we finally obtain 

limsuplimsupA' fc+fc _(^e) < W({Si:+ k ~}) (5.28) 
where, for any collection {M xy } of local functions on O x O, we defined 

W({M xy }) = sup I dx I dy{Q x ®Q y {M xy )-h{Q x \P)\ 

{Q x }CM T JT d JT d L } 

To remove the cutoffs we can argue as in the proof of Theorem 5.2 in [20], then we just 
sketch the argument. Let S a ^, = mm.{S aa i] k+} and define S xy and X k+ (£, e) accordingly. 
Using the boundness of S ^, ~ and that h(-\P) has compact level sets, one can prove that 
lim W({S x + h -}) < W({S X +}). Then, since X k+ {£,e) < X k+k _{£,e), from (5.28) we get 

fc_ ^oo 

lim sup lim sup X k+ (£,e) < W({S%+}) (5.29) 

£->oo ej.0 
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Now we are left with the upper cutoff. Recalling definition (5.24), by Holder inequality, 
for any p, q > 1 such that p~ l + q~ l = 1, and any k + G IN. 



X(£,e) < E ^ D ^ 1 log iE®^ exp 



+ 



p 



log-ZE 0CTjP<T exp 



? £<! El^!((^')rfc)J( w -^) 



(5.30) 



To bound the first term in the r.h.s. of (5.30) we can use (5.29) with S + , replaced by 



pS^, . To bound the second term we recall first that there are C > and a bounded subset 
B of IR d such that S aa >(uj u uj 2 ) < C(l + N B (wi) + N b (uj 2 )) for any a, a', so that 



(S^-^)(wi,W2) < [C^ + NbM + NbM) -k+] + 

<\[C + 2CN B - fc+] + + i [C + 2CiV B ( W2 ) - fc+] + 



Then, setting 



we have 



Yi,k+ M = exp 



/ rfr [C + 2CW B (r r w) - fc+] 



iE®^ exp 



IE 1 



< IE 1 



Y e,k+ M 



so that, from (5.30), 



f 



lim sup lim sup X(£,e)< lim W({p5j+})+ lim lim— — logiE' 



(5.3i) 



The second term in the r.h.s. of (5.31) is zero, see [20, Thm. 5.2]. Then, since the first 
term in (5.31) is bounded by W({pS xy }) , in the limit p I 1 we finally get 



lim sup lim sup X(£,e) < W({S xy }) 



»oo ej.0 



(5.32) 
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But (recall the above definition of S xy ) W({S xy }) is equal to the r.h.s. of (5.17) with 
fj, = 0, which then follows from (5.23) and (5.32). The theorem is proved. □ 

Now we can give the missing proofs of Section 3. 

Proof of (3.21). From the basic entropy inequality, for any 5q > 0, 



jfif-/ (*) 



dx P0, p (t, x) 



- S^H^t) < 6q 1 i d log IE exp 



5 -f~ d / dxf3tt p (t,x) 
Jr d 



(5.33) 

As before we forget the dependence on t G [0, T] and we deal with a generic smooth map 
x i— > X(x) G K u such that the corresponding density p(x) = d\p(/3, A(x)) G -£T TO satisfies 
(5.9) (the compact sets K u and K w have been introduced just after Theorem 2.1). Let / 7 
be the density of the corresponding local Gibbs state. From (5.33) it is enough to prove 
that, for Sq small enough and uniformly in the choice of x i— > A(x) above, 



lim sup lim sup lim sup 5 Q 7 log IE' 1 exp 



5 7" d / rfx/30 p (a:) 



0, p = 3,4 (5.34) 



where ft p (x) is defined as in (3.20) but relative to the map x 1— > A(x). We analyze the 
cases p = 3, 4 separately: 

(p = 3). From the smoothness of x ^ X(x) and since \<p * il n (u;. j7 )| < ||<£>||oo) f° r an Y 
(p G C(T d ) and any u;. j7 , there is Ci > such that 

fi 3 (aO < Ci(l + i? ) + (0fc O 

))d )) (5-35) 



where 



q,q £w 

q^q 



From Theorem 5.3, (5.11) and (5.35) we get the following bound: 



(5.36) 



lim sup lim sup lim sup 5 1 ^ d log IE Jl exp 



if-, 



5 o7 " d / dx/3n 3 (x) 



I J^dx [SoPdQ* ((1 + Rn + (</> k o |G|) n ) (1 - u n )) (5.37) 



< lim sup lim sup sup < I dx 

fc^oo n^oo [ JT d 



-Kp Mx) {Q x ) -(3p{p,\(x)) ; {Q x } C M T : / cfep(Q x ) = l 
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For any Q x G M T let J v x [de) Q e be its ergodic decomposition. Recalling that Kp^\(-) is 
an affine functional, the r.h.s. of (5.37) becomes 

limsuplimsupsup <^ / dx v x (de) 5 /3CiQ e ((l + R n + (<Pk ° \G\) n ) (1 - u n ) ) 

- K pMx) (Q e ) - Pp(p, A(x))] ; {zaJ : ^ d rf;r ^ ^(de) p(Q e ) = 1 



(5.38) 

Since |(l + i? n + o |G|) n )(l — u n )\ < l + k + R n and, for any {v x } in the supremum above, 
fdx J v x (de)Q e (l + k + R n ) < 2 + k, we can pass to the limit n — > oo inside the supremum 
and apply the Dominated Convergence Theorem. Then we can drop the constraint so that 
(5.38) is bounded by 

limsup / ctesupj fu(de) \s /3C 1 (l + p(Q e ) + Q e ((f> k o\G\))(l - l K {p(Q e ))) 

fc^oo JT d v y J L 

(5.39) 

-K PMx) (Q e )] - (3p{(3,X{x)) 

By arguing as in the proof of Lemma 4.2 of [25] (see also the proof of Lemma 3.1 in Section 
4), since V is positive and superstable, there is C 2 > such that 

\G\(u) <C 2 pU(u) (5.40) 

where U(oj) is defined as in (5.4) with an appropriate choice (depending on \G\) of the 
function \- Then, from (5.3) it follows that, for any Q G -M T , 

Q(<p k o \G\) < C 2 (3$(Q) 

so that, setting 5 = 5 max{/3CiC 2 ; Ci}, (5.39) can be bounded by 

J dx sup J u(de) {(6p + 6Pp(Q e ) + 6p$(Q e ) - K PMx) {Q e ) -pp(p,X(x))) 

X (1 - Mp(Qe))) " (W,A(X)) + ^,A(x)(Qe))lx(p(Qe))} 

< / cte sup {(5/3(1 + (A(x) + <y)m) -5), X(x) +5) 

Jq-d m ^ 

- pp{p, X{x)) - I p(1 _ s) (X(x) + 5, m)) (1 - l K (m)) - Ip(X(x), m)1 K (m)} 

In the last inequality we used (5.5), (5.8) and introduced the functional Ip(X, m) defined in 
(3.27). Observe that, for any A G K u , m \— > /^(A, m) is non negative, convex, and strictly 
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positive for m £ K. Recall also that the pressure p(/3. A) is a continuous function of its 
variables. Then, by continuity, the superior is for S (i.e. 5 ) small enough. 

(j) = 4). Since x h- > D^ rj X(x) is bounded and \D^ V J * R n (u. j7 )\ < \\DJ\loo, recalling 
definitions (3.12) and (5.36), there is C 4 > such that 

/ dxVt 4 (x)<C 4 l dx \\G\ n {u xn ) - k] + = C 4 dx \\G\{uj xn ) - k] + 

JT d JT d JT d 

Then (5.34) for p = 4 is proved if, for Sq small enough, 

= (5.41) 

Arguing as in the proof of Theorem 5.4, (5.41) follows if, for any 5 small enough, 

lim sup lim sup 7 d log JE^ exp [T (s > k) (u)] - [ dx 0p(p, X(x)) = (5.42) 

fc^oo 7J.O JT d 

where P 1 is the Poisson process on T d with intensity r )~ d and 

T^\u) = Y,PK*)-^ E V ^ z - *')) + 5l~ d f dx [\G\(u x „) -k} + (5.43) 



1/7 f 

lim sup lim sup S 7 log JE' 1 exp 

fc^oo 7J.O 



5 f3C 4l - d / dx[\G\(u xn )-k] 

JT d 



z,z Eto 
z^z' 



Let £q be such that U(u>) = U(u> fl Di ) with U chosen as in (5.40). As in the proof of 
Theorem 5.4 divide T d into disjoint boxes B a of side (2 7 £) and assume X(x) constant on 
this partition for small 7's. Let B a ^ be as in (5.18). From (5.40) we can estimate: 



T d J dx[\G\{u xn )-k] + <^ E Vi^iz-z')) 

T d \\JB^ lo 



(5.44) 



z,z'eu\UB^,2£ 
2/2' 



Restricting to 5 < C 2 1 , neglecting the interaction between different boxes and using 
the independence properties of P 7 , from (5.44) we get 

IE P ~< exp [T^ (u)] < IE P ~< exp [T^] (5.45) 



where 



2^2' 



+ 



E n7- 1 (^-^))+5 7 - d y ^[|GiK 7 )-fc] + 



2,2'6wn(B CT \B (T ,2^ ) 

2/2' 



Wo 
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Now, expanding variables, we can rewrite: 

JE^ exp [T^ k \oj)] = zg>(/3, X(a))IE< a) exp Id [ dr [\G\{r r u) - k] + 



(5.46) 



where is the grand canonical measure on Dg with chemical potential A(cr) and inter- 

action energy 

5C 2 



where -D^ = Dg \ Dg_ 2 g , the Hamiltonian Hg(uj) was defined in (5.2), and Z { *\[3, X(a)) 
is the corresponding partition function. From (5.45) and (5.46) we get 



7 



log JE P ^ exp [T<**>(u/)] < 1 d \Dg\ £ l^l" 1 logZgjGS, A(a)) 



(5.47) 



Since F is superstable and positive, [23], there is a positive constant b such that, for 



any u E O and any bounded region B C iR , 



q,q'euCiB 
q+q 



B 



Then we can estimate 



^S(^A(a))<iE^exp 



(3(1 - SC 2 )bNl o 
[3\(a)N Di _ 2lo - [3Hg_ 2 g e '*° + \{a)N Dl io 



\ D i,to\ 



< 



exp 



4/36(1 - 5C 2 ) 



\Dg,g \ 



so that 



limsuplimsup7 d | J D £ | V \Dg\~ 1 log zff (p, X(a)) < [ dxpp{p,X{x)) (5.48) 
i^oo 7 io „ ir d 



On the other hand, from (5.40), 



((5,<7) 

iE MD <? exp 



a / rfr [|G|(r r u;)-fc]' 



< IE M ^ exp [6C 2 /3He(u)] (5.49) 
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and, again from superstability, the r.h.s. of (5.49) is finite for any £ if 5 is small enough. 
Then, by the Dominated Convergence Theorem, 



lim sup IE D e exp 



k—>oo 



8 I dr [\G\{T r oo) - k] + 



so that 



lim sup lim sup lim sup 7 d | .D^ | ^ \Dg\ 1 \ogIE fJ ' D t exp 6 / dr [\G\(T r uj) — k] + 

t^oo k^oo 7J.O L JD e _ e 



= 



(5.50) 



From (5.47), (5.48) and (5.50) we get (5.42) (in fact the l.h.s. of (5.42) cannot be negative). 
□ 

Proof of (3.25). Recalling definition (3.20), by Theorem 5.4 the l.h.s. of (3.25) can be 
bounded by 



lim sup sup 



dx / dy 5f3[(\- f3\V\\ 2 )(s,x)Q x (R n u n ) -0A\(s,x) 

T d JT d L 

x Q x (P(R n )u n ) - /3VA(s, x) ■ VJ{x - y)Q x (R n u n )Q y (R n ) 

- PAJ(x - y)Q x (P(R n )u n )Q y (R n ) - fl(s,x, p(s))Q x (u n )] 

- K[3 M s,x)(Qx) - pp(P,X(s,x))];{Q x } C M T : / dx p(Q x ) = 1 



(5.51) 



For any Q x G M. T let f u x (de) Q e be its ergodic decomposition. Just as argued after (5.38), 
we can pass to the limit through the supremum and apply the Dominated Convergence 
Theorem. Then, calling v x (dm) the distribution of p(Q e ) under v x (de), (5.51) can be 
bounded by 



sup < / I j| h> z {dm(z)) I dx 5/3 [0(s, x, m) — 0(s, x, p{s))\ l^(m(a;)) 
17 z er d Jrd L 



Ip(\{s,x),m{x)) ; {v x } : dx v x (dm)m < 00 



(5.52) 



from which (3.25) follows immediately. □ 
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